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Abstract 

We present detailed computations of the nondecaying terms (three dominant orders) 
of the free energy in a one cut matrix model with a hard edge a, in /3-ensembles, with 
any polynomial potential. (3 > is not restricted to the standard values (3 = 1 
(hermitian matrices), /3 = 1/2 (symmetric matrices), (3 = 2 (quaternionic self-dual 
matrices). This model allows to study the statistic of the maximum eigenvalue of 
random matrices. We compute the large deviation function to the left of the expected 
maximum. We specialize our results to the gaussian /3-ensembles and check them 
numerically. Our method is based on general results and procedures already developed 
in the literature to solve the Pastur equations (also called "loop equations"). It allows 
to compute the left tail of the analog of Tracy- Widom laws for any (3, including the 
constant term. 
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1 Introduction 



Extreme value statistics, in particular the statistics of the maximum of a set of random 
variables, play a very important role as they have a large variety of applications in 
various fields such as finance, hydrology (rainfall, flood...) and engineering. On the 
other hand, random matrix theory has raised a huge interest in both physics and 
mathematics communities. The distribution of the maximal eigenvalue of a random 
matrix is thus of great importance. 

For usual Gaussian ensembles of random matrices (Gaussian Orthogonal Ensemble, 
Gaussian Unitary Ensemble and Gaussian Symplectic Ensemble), the distribution of 
the maximal eigenvalue close to its mean value is known to converge to the Tracy- 
Widom distribution [25, 26] for large N (where the matrix is of size N x N). This 
describes small typical fluctuations around the mean value. We are interested here 
in the large deviation of the distribution, which describes large and rare events (well- 
spaced from the mean value). We actually consider as starting point a distribution 
on M. N that generalizes the distribution of the eigenvalues of a random matrix for a 
general potential V(x) (instead of a quadratic potential for Gaussian ensembles) and 
a general /3 (Dyson index). 

We consider the following joint distribution: 



1 N 

P(A l7 ...,A^) = — -|A(A 4 )r l[e-^ v ^ with A(A,) = J] (A, - A,) 

^ ' i=l l<i<j<N 

(1-1) 

Z(oo) = J RN dAi . . . dAyv |A(Aj)| 2/3 Yli=i e is the (unconstrained) partition 

function. 



We define: 



Z(a)= [ dA 1 ...dA^|A(A J )| 2 ^ TT e -^<*> (1-2) 

J]-oo,a] N 7=1 

Zia) is the partition function for the weight given by Eqn. 1-1, multiplied by YLiLi ©( a— 
Aj). The probability that the maximal eigenvalue A max is smaller or equal to a (i.e. the 
probability that all eigenvalues A« are smaller or equal to a) is: 

V(a) = [ d\ 1 ... d\ N P{\ 1 , ...\ N ) = J^- (1-3) 

J]-oo,a] N Z y°°) 

The probability that the largest eigenvalue lies between a and a + da is: 

d a Vda = ^- (1-4) 
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For (3 = 1/2, 1, 2, the probability distribution on M N (the space of eigenvalues) given in 
Eqn. 1-1 derives from an invariant probability distribution on a set of matrices, whereas 
no such derivation is known for other values of (3. We still call Z(a) a "matrix model", 
and talk about " /3-ensembles" . This is only a convention of language to refer to a 
statistical weight of the form Eqn. 1-2. In the case of a quadratic potential V(x) = 
one recovers the usual Gaussian ensembles: Gaussian Orthogonal Ensemble ((3 = 1/2), 
Gaussian Unitary Ensemble ((3 = 1) and Gaussian Symplectic Ensemble ((3 = 2). And 
in general, one defines the Gaussian (3 Ensemble. Note that (3 differs from the usual 
Dyson index (3 by a factor 2: (3 = 2(3. The normalization we choose in the definition 
of G/3E is also the one used in [24]. Let us mention Dimitriu and Edelman [13] have 
found a tridiagonal matrix model, whose distribution of eigenvalues reproduces the 
G(3E ensemble. As far as our methods are concerned, these sparse matrix models are 
very different in nature from those we consider. 

We want to compute the large A" expansion of Z(a). For a one cut solution, we 
assume 1 that lnZ(a) has an expansion in powers of 1/N. It takes the form: 




coefficients F 9,k are independent of A" and (3. They may depend on a, on t, and on the 
potential V(x). They are given by geometric quantities associated to a plane curve, 
which were introduced in [11]. In this case, the plane curve is related to the equilibrium 
density of eigenvalues. Dx,p is a normalization constant, it depends only on A" and (3, 
and not of a,t and V. We tackle the question of this constant in Section 3.8. 



For instance, limiting ourselves to the non-negative powers of A" (we will call them 
"nondecaying terms" or "stable" terms), we have: 

N 2 N 

lnZ(a) = D Ntf) + — /3F ' (a) + -(/3-l)F ' 1 (a) + F 1 ' (a) + (/3 + r 1 -2)F ' 2 (a) (1-6) 

F°'° is called the prepotential, F 0,1 is the entropy, F 1,0 is related to a Laplacian de- 
terminant, and F ' 2 is the Polyakov anomaly. All the higher F 9,k are also given by 
geometric quantities (their interpretation is less known), and there exists an algorithm 
to compute them recursively. 

For the statistics of the largest eigenvalue, F°'°(a) is the function of large deviation, 
and higher F 9,k are its subleading corrections. This expansion of F(a) as a series in 
1/N for large A" is actually valid only on the left of the mean value. On the right 

1 The proof of this assumption under some hypothesis is a work in progress [9] 
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of the mean value, In (Z(a)/Z(oo)) is expected to be exponentially small, i.e. much 
smaller than any power of 1/N. Note that, for the Gaussian case (quadratic potential 
V(x) = Tp), the leading term F°'°(a) of the left-tail large deviation has already been 
computed a few years ago [14, 15] using a Coulomb gas technique [23]. Here, we explain 
how to compute recursively all the F 9,k for a general potential and a general /3. We 
also give an explicit expression for the nondecaying terms F°'°, F ' 1 , F 1 ' and F 0,2 
for the Gaussian case (quadratic potential) and compare our results with Monte Carlo 
numerical simulations. 

Outline of the article 

We first present (Section 2) the topological expansion of matrix models and the method 
of loop equations. We describe in details the one-cut solution of the loop equations in 
Section 3. We specialize it for the Gaussian-/? ensembles in Section 4, with comparison 
to numerics in Section 4.2. We discuss in general the scaling limit at the edge of the 
spectrum, and make the link with Tracy- Widom law in Section 5. The appendices give 
the proof of the results used in the text to settle the "constant problem". 
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2 Correlators and loop equations 
2.1 Definitions 

We shall consider the expectation value of the resolvent 

and more generally, the non-connected correlation of n resolvents: 

n N 



W n (x u ... t x n )= (f[£—L-) (2-2) 

j=l i J= l 3 % . 



Their cumulants (the connected correlation functions) are called for short "correlators": 

I n N \ 

W n ( Xl ,...,x n )= (IIE^A - / ( 2 -3) 



\j=l i j= l o 



c 



For instance: W 2 (xi, x 2 ) = W^Xi, ^2) + ^1(^1)1^1(^2) 
Let us write the potential: 



V(x) =t + J2~ x3 ( 2 - 4 ) 



and introduce the loop insertion operator: 



Applying successively this operator to lnZ(a) gives the n-point correlation functions: 

w) ww)'--ww) lnZ[a) (2 ' 6) 

The correlators satisfy loop equations. These are named Schwinger-Dyson equations 
(in physics literature), or Pastur equations (in maths literature). They follow from the 
invariance of an integral under a change of variable. Alternatively, they can be obtained 
by integration by parts. One must just pay attention to boundary terms at A, = a. 
We derive them by two equivalent methods, giving complementary information. 
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2.2 Loop equation - 1st method 

Consider an infinitesimal change of variables Aj — > A j + e f^f^ . This change of variable 
preserve the range of integration ] — oo, a]. We have, to first order in e: 



A 



A 



i=l i=l > 

|A(A t )| 2/3 |A(Ai)| 2/3 
|A(A 4 )| 2 ^ 



x — a 



(x - \y 
1 + ^ E 

l<j^j<AT v 

1 + J] 



+ 0(6 2 ) 



x — a 



Aj x — Aj / Aj — Aj 
+ 0(e 2 ) 



0(e 2 



x — a 



l<j^i<N 



n 



N/3 

e * 



V(A0 



n 



(x - Aj)(x - Aj) 



1 ^^ nAQCV-a) | Q(e2 , 



JV/3 



J 2 W 



t=l 



x — A,- 



(2-7) 
(2-8) 



Writing that <5 e In Z(a) = 0, we get to order 1 in e: 



t 2-^ x - A j 

i=l 

N 

Ex — a y—^ 
, (x-A,) 2+/ ^ ^ 
i=i v ' 



x — a 



l<j^i<N 



X — 



1-/3 



(x - Aj)(x - Aj) 

N „ AT 



= (x - a) ((/3 - 1) W{(x) + /W x (x) 2 + (3W 2 (x, x)) 

where we have used the notations of Eqn. 2-3 above. 
The left hand side can also be written: 



(2-9) 



N 

E 

i=l 



V'(Xi) (A, - a) 

X - A,; 



N 



(x - a) 



t=i 

N 



V'(Xi) 

X - A,; 



iV 



t=l 



X 



V"(A,) - V"(x) + V'(x) 



N 



i=l 



x - Xi 



i=l 



A' 



(x - a)V"(x)iyi(x) - (x - a) ( 



nx) - nAj ) 

x - Xi 



N 



Em) 



(2-10) 



t=l 



Let us define: 



N 



X - A,; 



which is a polynomial of x, of degree cleg V — 1, and 



N 



£m 



(2-12) 



i=l 



The loop equation can thus be written: 



PW^x) 2 + PW 2 (x, x) + {P~ l)W[{x) 



N/3 



V'(x)W 1 (x)-P 1 (x) 



x — a 



(2-13) 



2.3 Loop equation - 2nd method 

We can also perform a infinitesimal change of variable changing the bound of integra- 
tion: 



A; ->■ A; + e 



+ 0(e 2 



x - Aj 

We have to first order in e: 

N N / N \ 



(2-14) 



i=l 



A(A) 2/3 -> A(A) 2/3 
A(A) 2/3 



^ la; — Aj x — A, J / 



1 



0(6 2 



iV 



l<i^j<N 

(2-15) 



v v 



n 

i=l 



i=l 



i=i 



i x — A,- 



(2-16) 



Under this change of variable, the integral changes by e °_ — + 0(e 2 ). Collecting 
the first order in e, we obtain: 

PW 2 (x, x) + /W(x)) 2 + {(3 - l)W[{x) - ^ (V'(x)W 1 (x) - P x (x)) - 



x — a 
(2-17) 



This is again Eqn. 2-13, with the alternative information that ^-c = —d a F(a 



2.4 Higher order loop equations 

By application of C^j ^^y, • • • , C^j dv ^ n _ l} to Eqn. 2-17, one can obtain loop 
equations for the n-point correlators. We have to introduce: 



P n (x , Xi, • • • , x n 



Tr 



V'jx) - V'{M) 
x - M 



n— 1 -. 



i=i 



x 7 - M 



(2-18) 



Let us write / = {x±, . . . , x n -i} for the spectator variables. The result is: 

dW (x I) 

£ PW w+1 (x, J)W n ^ A {x, I\J) + 3W n+1 (x, x, I) + (/3- 1) 

jci ax 

= ^ ( W „(,, /) - P„(x ; /)) + E T" ( ^-^Ml-wwm 

* (2 . 19) 

x — a 

2.5 Recovering F 

The easiest way to find F is to integrate its derivative with respect to the independent 
parameters (a,t,tj). These derivatives can be expressed in terms of the correlators. 

• Wrt the hard edge a. By comparison of the loop equations obtained by the two 
change of variables, we have seen that: 

d a F = -— (TrV'(M)) = — Res dx V'(x)W 1 (x) (2-20) 

t \ it x— >oo 

This could also be obtained by comparing the 1/x term (when x — > oo) in 
Eqn. 2-17. Another method uses the fact that Pi is a polynomial, thus regu- 
lar at x = a. So, we have: 



d a F = lim(x — a) 



8W 2 (x, x) + BW^xf + (B- l)W[(x) - —V'{x)W 1 (x) 



(2-21) 



We will see later another method based on a residue formula. 
Wrt the times t and tj. We have from the matrix integral: 



NB 
1 

NB 



dtF = -f(V(M) 



t 



Res dxV(x)W 1 (x) (2-22) 



x— >oo 



V - -f <? ) 



t x— >oo j 

2.6 Topological and h expansion 
2.6.1 Definition 

We introduce two parameters: 



NB x j 

Resdx— W 1 (x) (2-23) 



t ' N \ 8 
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We shall consider the large N expansion of the free energy F(a) and of the correlators 
W n (xi, . . . , x n ), in the form: 

F = v 2 ~ 29hk F9 ' k ( 2 - 25 ) 

g,k>0 

W n {x 1} ...,x n ) = f3- n ' 2 v 2 ~ 2g - n h k Wi> k (x 1} ...,x n ) (2-26) 

g,k>0 

The prefactor in Eqn. 2-26 was chosen such that: 





dV{x n 



W£ k ( Xl , ...,x n ) = W&fa, . . . ,x n ,x n+1 ) (2-27) 



For the /3-matrix ensembles, h is of order 1/N. Those definition are such that F 9,k and 
the functions W%' are independent of N and f3. They depend on a, on t, and on the 
potential V(x). For instance, the free energy is: 

F(a) = ^-/3F^ + -((3-l)F^ l + F 1 '° + (/3 + (3~ 1 -2)F^ 
t 2 t 

+ ^ ((1 - fi-^F 1 ' 1 + ((3-3 + 3/T 1 - /T 2 )F°' 3 ) 

+^ (p-^F 2 ' + (1 - 2/T 1 + p- 2 )FW + ((3-4 + Q/3- 1 - 4/T 2 + /T 3 )F ' 4 ) 
+o(l/iV 2 ) (2-28) 

2.6.2 Expansion of the loop equations 

All the same, we expand: 

oo 

P n (x,I) = (3- n ' 2 v 2 ~ 29 h k P%' k (x,I) 

g,k>0 

To the leading order, the loop equation (Eqn. 2-17) reads: 



(W?'°(x)Y = V'(x)W?'°(x) - P^(x) - (2-29) 

Ob Qi 

This relation is sometimes called master loop equation. And, for higher g, k, we 
have: 

(V'(x) - 2W?'°(x)) Wt k (x) = Wr 1>fc (x, x) + k -\x) 

ax 

i 



+ J2 wf' k '(x)wr 9 '' k ~ k '(x) 

0<g'<g,0<k'<k 

+Pf(x) + — (2-30) 
x — a 



And, for n-point correlators, the expansion of Eqn. 2-19 is: 



(V'(x)-2W^(x))W^(x,I) 
W^ 1 1 ' k (x,x,I) + ^-W 9 ' k - 1 (x,J) 



+ J2 W°/ +1 (x,J)W^j\ k - k (xJ\J) 

o<g'<g,o<k'<k,jci 

w^iMx^-^wtMi) 




+P^\ X J) + -^±1 (2-31) 

x — a 

In those expressions, Yl' m eans that we exclude of the sum the terms where W®' 
appear (we put them by hand in the LHS). 

The general solution of those equations was found in [11, 12]. All W 9,0 are "ge- 
ometric covariants", and all F 9,0 are "geometric invariants", associated to the plane 
curve C: 

(x,y)eCCC 2 & y = YM^ W W( x) (2-32) 

The function y defined on £ is closely related to the density of eigenvalues in the ther- 
modynamic limit (N — > oo). The triplet (x,y,C) is often called the spectral curve. 
A deep property of stable F 9,0 is their invariance under symplectic transformations, i.e. 
all transformations (x,y) —> (xi,yi) such that |dxAdy| = \dx\ Adyi\. It is conjectured 
that stable F 9 ' k for k ^ are also invariant in this sense, but there is no proof at the 
moment. So, we often call F 9,0 the "symplectic invariants of the curve", but we use 
the name "geometric quantities associated to the curve" to refer to general F 9,k . 

In practice, the first step is to compute the unstable correlators, that is the corre- 
lators W 9,k such that 2 — 2g — k — n > 0. The most important one is the resolvent to 
leading order, W 1 ' . The other unstable correlators are W^' and W"'\ Then, all the 
remaining correlators follow. In a similar way 2 , there are general formulas for unstable 
free energies (F°'°, F 0,1 , F 0,2 , F 1,0 ) and a uniform algorithm to compute all the other 

Notice that unstable terms come with a nonnegative power of N in the topological 
expansion, whereas stable terms are coefficients in the o(l). In a sense, this o(l) is 
easier to compute than the nondecaying terms. 

Let us also mention that all W 9,k (and Wq ,k = F 9,k ) have a combinatorial interpre- 
tation, as counting (possibly non-orientable) ribbon graphs of genus g, with k Mobius 
strips, and n marked vertices [12]. 



2 We may identify the free energies F 9 ' k to "0-point correlators, W$ k " . 
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.3 First few loop equations 
• For n = 1: 

2y(x)W?' 1 (x) = -^-W?'°(x) + P?' 1 (x) + -^— (2-33) 
ax x — a 



2y(x)Wl'°(x) = w2> (x,x) + P}'°(x) + -?— (2-34) 

x — a 



2y{x)W° 1 '\x) = A.W?>\x)+ «' 1 (x)) 2 + P 1 °' 2 (x) + ^- (2-35) 
and so on. 
For n = 2: 

2y(x)W°'\x, x') = -^W 2 °' (x, x') + W%'\x) ( 2W°'°(x, x') + 



dx \ (x — x') 2 

d fx' -a W°'\x')\ i p0jl ^ , c£'V 



dx' \ x — a x — x' J x — a 



2y(x)lF 2 1 ' (x,x') = W3 -°(x,x,^) + JOs) (W^x') + ) 



d ( x> - a Wl'\x>)\ , p i,o, n + c 2 '°(x') 



dx' \ x — ax — x'y x — a 

2y(x)lF 2 ' 2 (x,x') = ^PF 2 °' 1 (x,x / ) + 21F 1 ' 1 (x)PF 2 ' 1 (x,x / ) 



+2<' 2 (i) hy 2 '°(x,x' 



1 



(x — x') 2 



x'-a< 2 (x'A , p o,2, „ , 4 2 (x') 



dx \x — ax — x j x — a 

and so on. 

For n — 3: 

o \n/0,0/ \ \TJ/0,0/ \ , W / 2'°( a: 5 ^2) . W%'°(x,Zi) 

2y(x)W 3 ' (x,xi,x 2 ) = 2W 2 ' (x,x 1 )W 2 ' (x,x 2 ) + — - + 



dxi x — Xi dx2 x — X2 

-P 3 ' (x;x 1) x 2 ) + C °' 0(Xl ' X2) 
x — a 



x — X1) 2 (x — x 2 ) 2 
W 2 °' (xi^ 2 ) 



and so on. 
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2.7 Density of eigenvalues 

The expected density of eigenvalues is denned as 

N 



= ( - A,) ) (2-36) 



i=l 



The linear form / — > f pf, defined on an appropriate space of test functions, has a 
large N expansion (in the weak topology) of the type 3 



*>*^E(¥)"" (^) ™ ^ 

The 1-point correlator, W\(x), is the Stieltjes transform of the density This means 
that, if T> contains the eigenvalue support, 

Wl(x) = I *U&1 (2-38) 
Jv x-x' 

Conversely, one can see the density as the discontinuity of Wi(x): 

2iiTg(x) = W 1 (x- i0 + ) -W 1 {x + i0 + ) (2-39) 
In the thermodynamic limit (N — > oo): 

- q(x) ~ e°'°(x) = — (W?'°(x + *0+) - W?'°(x - t0 + )) (2-40) 

N N^oo 2Z7T V ' 



The origin of T> will become more precise when we present the analytical structure of 
W®' derived from the master loop equation. 

When x — > oo, we know that W^ ~ t/x, and this implies the normalization: 

dxg°'°{x) = t (2-41) 

v 

3 One-cut solution to the loop equations 

The spectral curve is defined by the master loop equation, that we rewrite: 



<'V) = ^ - J^f- - " ^ (3-1) 



3 In general, this is not true when p(x) is considered as a function. p(x) has a pointwise p°, but 
exhibits beyond leading order fast oscillations, even in the one cut case. We thank a referee for 
pointing an erroneous statement in an earlier version of the text. 
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If eigenvalues accumulate in the large N limit on one segment, of the form [b,a[, we 
know that Wi°'°\x) must have only one cut [b, a] in the complex plane, and thus 
W 1 ' (x) must be of the form 



r0fi ,^ V'(x) M{x) x 



Wr(x) = -p--p^— (3-2) 
2 2 V x — a 

where M(x) is a polynomial to determine. 

In principle, knowing that Wi'°(x) has one cut and matching Eqn. 3-2 with the 
large x behavior W^ ,0 (x) ~ t/x, allow to compute explicitly the endpoint b, as well as 
the polynomials M(x) and P 1 0,0 (x). The easiest way to do this computation is to use 
Zhukovsky map. 

The purpose of this article is to illustrate our methods, so we assume that: 

• the support of eigenvalues is of the form [b,a[, in particular it is connected 

• [b,a[C R 

Then, one can prove that W%' k (x±, . . . ,x n ) has one cut on [b,a[ in each variable for 
Xi G [b, a[, and that it is holomorphic in each variable outside this cut. 

We present detailed computations of the unstable correlators (W^' , W®' , W®' ) 
and the unstable free energies (F°'°, F ' 1 , F ' 2 , F 1 ' ), and the algorithm to compute 
the other W-[ ,k and F 9,k . General formulas for W^ and the procedure to compute all 
W k,g when the support of eigenvalues is a finite set of paths in the complex plane can 
be found in the literature [10, 11]. 

3.1 Spectral curve, W^ and density 

The Zhukovsky map is defined as: 

x(z) = a + 7(2 + 1/z) (3-3) 
with a = (a + b)/2 and 7 = (a — 6)/4. We also have: 



7(2; -1/z) = y/(x- a)(x - b), J 



x — b , 1 + z 



x — a 1 

r0,0, 



(3-4) 



It is then clear that W 1 ' (x(z)) and p[x{z)) are rational function of z (for this reason, 
z is called a uniformizing variable). 

Consider any potential V of finite degree d: 

d t 

V(x) =t +>4x J (3-5) 
1 

3=1 J 
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We have: 



V'(x) ~ t d x d ^ (3-6) 

x— ¥00 



We assume that V is such that the matrix integral converges. Let us decompose 

d-l 



V'(x(z)) =u + J2Mz k + 



(3-7) 



k=l 

With the Zhukovsky variable, W\(z) = Wi'°(x(z)) is the solution to the problem: 
(z) Wi is a meromorphic function defined on C. 
(ii) W\ is meromorphic outside the unit disk, and : 

Wi(z) ~ t/x(z) 

z—too 

(in) w\ is regular at z — — 1 (i.e x = b), and has atmost a simple pole at z — 1 (i.e 
x = a). 

(if) twi satisfies the following equation (deduced from the loop equation for a 1-cut 
solution): 

Vz e C, wi(x(z)) + = V'(x(^)) 

The solution is unique, and it is easy to find the answer in the form: 




The behavior when z — > oo determines b through: 

uo = -r , 7(r + «i) = t 
Then, the density in the thermodynamic limit p(x) = p 0,0 (x) is: 



(3-8) 



(3-9) 



p(x(z)) 



2in 



{wxil/z) - w 1 (z)) 



l + z 



d-l 



2in \ 1 — z 



+ ^2u k (z k - z 



(3-10) 



fc=i 
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3.1.1 More on the density 

It is sometimes convenient to use the variable 9 such that z = e ld , that is x = a + 
27 cos 9. We have: 

p(x(e ie )) = ^-(r cotan(fl/2) + 2^u k sin(Jfe0) J (3-11) 



k=l 



Eqn. 3-7 is in fact a decomposition of V'(x) on the basis of Chebyshev polynomials of 
the first kind T n (cos9) = cos(n#). 

V\x) = u + 2u k T k (cos 9) (3-12) 

k>l 

Likewise, it is useful to decompose V on the basis (T n ) n6N : 

V(x) = v + J2 2v k T fc (cos 9) (3-13) 
fc>i 

Then, V' is naturally decomposed on the basis of Chebyshev polynomials of the second 
kindf/„(cos£) = ^^. 

V'(x) = -Y j kv k U k (cos9) (3-14) 

1 k>l 

We recall that T n are orthogonal polynomials for some scalar product (-|-)r: 

1 f n 1 

(T m \T n ) = - d9T m (cos9)T n (cos9) = -5 nm ifra^O (3-15) 
7T Jo 2 

— 5o, m if n = (3-16) 

Similarly, U n are orthogonal polynomials for a related scalar product ( • | ■ ) : 

{U m \U n ) u = - T d9sin\9)U n (9)U m (0) = ^„, m if n ^ (3-17) 
71 Jo * 

= 5 ,m ifn = (3-18) 
The expansion of Eqn. 3-14 is related to Eqn. 3-7 through: 

\/k>2 j(u k -i - u k +i) = (Tk-i - T k+1 | y^nv n U n ^ 

n>l 

= ^2 sin 2 6> ?7 fe _i | ^nv n U n _^j 



n>l 



n>l 

kv k 
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u 



One can check that this formula holds for k = 1 as well. 



Vfc > 1, 



7(/Ufc_i - itfc+i) = kv k 



(3-19) 



So, we have a compact expression for the density dx p(x): 




(3-20) 



The eigenvalue support T> = [b, a[ is mapped to the upper unit circle from — 1 to 1 
(we call it — C + ), and to [ir, 0[ in the variable 9. The —1 sign in Eqn. 3-20 is consistent 
with this reverse orientation in the variable 8, since the total weight of p is t. 

3.1.2 The moment polynomial 

The density can also be written: 



From Eqn. 3-1, one knows that p{x) behaves as yx — b when x — > b. So, M must 
be a polynomial, and it is called the moment. Since the support of eigenvalues stays 
connected, b e] — oo, a[ is determined as the unique zero of the polynomial in bracket 
such that p(x) remains nonnegative for x G [b, a[. Besides, the degree of M is d — 1, so 
it admits d — 1 zeroes X\, . . . ,Xa-i- For generic V (i.e a non critical potential), those 
zeroes are distinct from a and b. They must lie outside ]a,b[, and are either real, or 
come in pairs of complex conjugate values. Let us call Sj, the point such that x(sj) = Xj 
and \sj\ > 1 (we also have x(l/sj) = Xj). We have: 




(3-21) 



d-1 



M{x) =t d Y[(x-x 3 ) 



(3-22) 



Also, we may write in a factorized form: 



tdTf' 



d-1 




p(x(e»)) 



cotan(#/2) ] 



(3-23) 
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For p e N, we may also introduce the following quantities: 

m (p) = Y t k?u k , =5](-lfFu fc (3-24) 

k>l fc>l 

They allow to relate the values of the p-th derivative of M(x) at points a and b, to the 
coefficients u k - By expanding Eqn. 3-25 near z — > 1 (x = a), or z — > — 1 (x = b): 



r 



M(a) = r ( M(b) = 4m^ 

M'(a) = 771^/7 S M'(6) = -(m[ 1} +2m[ 3) )/37 



(i) 

6 



3.1.3 The spectral curve 

The function y defining the spectral curve is given by: 

y{z) = inp(z) = -(w 1 (l/ z)-w(z)) 

1 M{x{z)) 



2z-l 

Y,u k {z k - z- k )\ (3-25) 



1 / z + 1 
- r 

2 z-1 



k>l 



3.2 The prepotential, F° 
3.2.1 General formula for F 

The prepotential may be found by integrating Eqn. 3-29, which might be useful for 
simple examples or numerics, but do give a hint for a general formula. For this purpose, 
one has from homogeneity (invariance of F with respect to (t, tj) y (At, At,) for A > 0): 

= I Res dxV(x)W°'°(x) + J<9 t F°'° (3-26) 
2 x— s>oo 2 

We have introduced the chemical potential \i = d t F ' . We state [19] that it can be 

computed as: 

u = It In 7 + Res V(x) (3-27) 

x-s>oo ^(xj 

With our notations, we find fi = 2thi7 — vq. 

To compute the first term in F°'°, we can move the contour from oo to C surrouding 
the cut. We have: 

Res dxV(x)W°'°(x) 

x—>oo 

= Res — lv + v k{z h + z~ k ) J ( r(l + 1/z) + ^ fc2 ;- fc J 
fc>i 
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Eventually, we collect the terms 



k>l 



(3-28) 

3.2.2 Computation of d a F° -° 

After Eqn. 2-21, d a F 0,0 appears as the coefficient of the simple pole l/(x — a) in: 

[«'°(x)) 2 - V'(x)W?'°(x) - P?'°(x)] 

Since is regular at x — a, and W^ 0,0 behaves as (x — a) -1 / 2 , we have when x — > a 
for z — > 1): 



«'°(*)) 2 



d a F 



o.o 



x — a 



d a F°'° 1 
7 (^-l) 2 



With Eqn. 3-8, we obtain: 



Q aF ^ = lr 2 



3.3 Finite correction to the 1-point correlator, W^ 1 

The loop equation (Eqn. 2-33) is: 



M{x] 



x — b 



Wi'\x) = ^W^(x) + P?>\x) + 



0,0/ 



30,1/ 



x — a 



,0,1 



Hence: 



W°'\x) 



dx 



M'(x) 



x — a 



Q(x) 



4(x-a) 4(x-6) 2M(x) M(x) y/(x - a)(x - b) 



(3-29) 



(3-30) 



where Q(x) is a polynomial. In the Zhukovsky variable z, this expression is clearly a 
rational function of z. We define: 



w?' 1 ^) = W^(x(z))x'(z) 



0,1/ 



(3-31) 



We see that a;°' 1 (^) has simple poles at z — 1 (with residue 1/2), at z = —1 (with 
residue —1/2), and at z — 1/sj, z = Sj and z — 0. Since W"' 1 ^) is holomorphic outside 
[6, a[, must be meromorphic outside the unit disk. Therefore, the polynomial 

Q(x) is determined such that the poles at Sj cancel. This implies that the poles at 
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1/sj have residue —1. Besides, (x) G 0(l/x 2 ) when x — > oo imply that u\' l (z) G 
0(l/,2 2 ) when z — > oo. This determines the residue of the pole at z — 0. Eventually 
we get: 



(3-32) 




3.4 The entropy, F° 1 
3.4.1 General formula for F 01 

The first correction F 0,1 was computed by Dyson himself [18]. In our notations, it is 
the opposite of the entropy: 



^0,1 



dx p(x) ln(p(x)/t) 



Let us insert the expression the density (Eqn. 3-20), and its factorized form 
(Eqn. 3-23): 



M,l 



— ( t - kv k T k [c 



k>l 



d-l 



d-l 



2nt 



cot an 



(0/2) ] 



3=1 



t In 



d-l 



2nt 



d-l 



(3-33) 



where: 



1 



n d9 



7T 



it - ^2kv k T k 

\ k>l 



[cos 9) ln|tan(0/2)| 



JU) = [ — \t-J^kv k T k (cosd) ) \n\e l6 -s\ 

To compute X, we use the classical integral J^ln |tan(6*/2) | = and do an integration 
by part: 

sm(k8) 



1 



n d9 
— 2_^v k 



7r * — ' sin 6 

k>l 



k odd 



We compute J(s) for \s\ > 1: 



J{s) 



/ — \ t-Y, kv k T k(^^) 

\ k>\ ) 

thx\s\+ I — It - ^2kv k T k (cos6)\ 

\ k>l J 

/ oo 

■ l-E Re ( 



In \s\ + X - ln(l - e i9 /s) + X - ln(l - e~ id ft 



, m ^ r m (cosg) | Im ^-m^ sin M) 



m=l 



?7? 



??? 



1 / °° 1 \ 

k>l k>l \m=l J 



A m ^ is a real coefficient which vanish when m and k are of the same parity, and else 
is equal to: 



1 1 

71 \m + k m — k 
We compute J7"(l/s) for |s| > 1 in the same way: 



d9 

7T 





77 



t - y^fct> fc T fc (cos 

fc>l 

I t - y^fcf fe r fc (cos6i) J 
\ fe>l / 

00 

E -Re( S " 



Ii n (i- e -7 s ) + Ii n (i- e ^/i) 



T m (cOS#) _ Im ( s -r^ Sm M 



s.m=l 



^^Re(s V^^Im(f> m ,fc^) 

fc>l fc>l \m=l / 



In the sum J{s) + J7"(l/s), the term containing an imaginary part cancel. Moreover, 
we know that if Sj is a zero of p, then ~sj is also a zero of p. So, . s~ fe is already real, 
and the absolute value in the log only remove the sign of Sj. Collecting the terms: 




(3-34) 
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3.4.2 Computation of daF ' 1 

One may also compute d a F 0,1 , which requires less algebra. We have when x — > a (or 
z — > 1) by expansion of Eqn. 2-21: 

daF ' 1 



[(2<»°(x) - V{xj) W?'\x) + d x <'°(x)] 



x — a 
d a F^ 



7 (z-l) 2 

We may extract d a F 0,1 by computing a residue. Since ~ 27 (2 — 1) when z — > 1, 
we write: 



- Res dzx'fz) 
2 w 



Res dz y(z)ui 1 ' (z 



(2wr( Z )-V'(x( Z )))w Az) + ^^ 

0,1/ 



(3-35) 



With u ' 1 given by Eqn. 3-32, and y(z) = iirp(z) given by Eqn. 3-25, we obtain: 



d^ ' 1 = Res ( — ^ + r - + 0(z-l 



-+1 U-l 2 



>)(BH-g^<H 



Hence: 



j'=i J 



(3-36) 



3.5 Leading order of the 2-point correlator, W 2 '° 



By application of the loop insertion operator to Eqn. 3-8, we find that w 2 (zi,z 2 ) = 
W 2 '° (x(zi) , x(z 2 )) is the unique solution to the following problem: 

(i) W2(zi, z 2 ) is a meromorphic function in z\ and z 2 , which is symmetric in its two 
variables. 

(ii) w 2 (zi, Z2) is meromorphic for z\ and z 2 outside the unit disk, and behaves as \ 
when z\ — > 00. 

(in) For z 2 7^ ±1, ^2(^1; ^2) has atmost simple poles when z\ — > ±1 (i.e when X\ — a 
or x\ = b). 

(iv) w 2 satisfies the following equation (deduced from the loop equation for a 1-cut 
solution): 

1 



Vzi e C, Vz 2 e C\ {±1}, w 2 (z 1 ,z 2 ) + w 2 (l/z 1 ,z 2 ) 



(x(z 1 ) - x(z 2 )) 2 
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This problem is the same with or without hard edges, and its solution is universal (it 
only depends on the endpoints a and b, through the Zhukovsky map), it is given by: 

1 1 1 



w 2 (zi,z 2 ) 



x'(z 1 )x'(z 2 ) (z 1 - z 2 ) 2 (x(zi) - x(z 2 )y 



(3-37) 



or equivalent ly: 



W 2 [Zt,Z 2 



x , (z 1 )x , (z 2 ) {z\Z 2 - ly 



(3-38) 



3.6 Unstable W^ g 
3.6.1 Analytical structure 

Initially, W 9,k (x(zi), . . . , x(z n )) was defined and meromorphic for \zi\ > 1 and Zi ^ 
±1. We shall see that it can be analytically continued (this is true under the 1-cut 
assumption) as a meromorphic function on the whole complex plane, 1.6. £LS db rational 
function of its variables. We call this function w 9 l ' k (zi, . . . ,z n ). Let us introduce 

L0 9 n ' k ( Zl , ...,z n ) = x'{ Zl ) ■ ■ ■ x'{z n ) w 9 n '\ Zl , ...,z k ) (3-39) 

We shall now investigate the analytical properties of . We already know the 
important property: 

(ii) uj% k {z\, . . . , z n ) has no singularity for \Zi\ > 1 and Zi ^ ±1. 

Let us rewrite the loop equation by separating "singular" and "regular" terms. We 
recall the notation: 



and we introduce 



w°'°(x{z)) = 
u g n ~ 1 {\z,z,I) + 



V'(x(z)) 



y( z ) 



-\A 



+ U \j\+l( Z > J ) u9 n 
JQI,0<g'<g,0<k'<k 

+d z {uf-\zJ))- X -^-u^\zJ) 



9 - 9 '> k - k '[z,I\J) 



+ J>'(z))X 



x'(z) x(z)—a x'(zi) 

x(z) - x(Zi) 



S 9 /(z,I) 

K k (z,i) 



2y{z)x'{z) 



x(z) — a J 2y(z) 



(3-40) 
(3-41) 
(3-42) 



Zi£l 
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The loop equations may be rewritten: 



u% k (zJ)=S°> k (z,I)+'R% k (zJ) 



(3-43) 



S^' h has the same polar structure as cu^' fc , apart from the fact that it may have a 
pole in the physical sheet, at z = Sj. To find u^' k , we should add lZ% k to S n ' k , thus 
cancelling this pole. The remainder TZ^ k (z,I) has the following properties: 

• TZ^ k (z, I) is regular when z — > ±1. 

• 7Z^ k (z, I) may have poles only at z = Sj, z = 1/sj, and they are atmost simple. 

By a recursion of 2g — 2 + n + k using the loop equations, It is then possible to 
complete the description of the analytical structure of uj g n k : 

(i) u^' k (zi, . . . , z n ) is a rational function of all its variables (thus defined on C n ). 

(Hi) u^' k (zi, . . . , z n ) may have poles at z\ = {±1}, z% = 1/sj (1 < j < d— 1), z\ = 1/zi 
(i ^ I), and Zi = 0. 

(iv) uj% k (zi, . . . , Zn) G 0(1/ zf) when Zi -> oo. 

It is clear that the loop equations, the analytical properties (i) — (iv), and the data 
of u>i'°, Ui' 1 , (J^° determines uniquely the family of rational functions oo!^ k . 

(v) We know an extra property of the id n ,k we are looking for: they are symmetric 
functions of all their variables. Although not obvious, the fact that this unique 
solution is symmetric can be proven from the loop equations. 

(vi) It is also possible to prove by recursion that the pole at z = is atmost simple 
(we assume d > 2). It comes from the leading term of P% ,k (x(z), I) (in TZ^ k (z, I)) 
when z — > (x — > oo), all the other terms in S g,k (z, I) and TZ^ k (z, I) are regular 
when z — > 0. 

3.6.2 Algorithm 

Let us recall that x(z) = x(l/z) and y(z) = —y(l/z). Since P^' k (x(z), I) is a polyno- 
mial in x(z), we know that: 



Since this rational function in z has only simple poles at z = Sj and z — 1/sj, it must 
be of the form: 



n k (i/z,i) 



z 2 n k ( z ,i) 



(3-44) 




(3-45) 
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To cancel the pole at z = Sj in S^ ,k , we must have: 



Resd£<S^(£,I) (3-46) 



Since 5^ ,fc (z, J) has only a simple pole at z = Sj, we have: 



n 9 ' k (I) = — { 3 ' (3-47) 
2y>{s j )x>{s j ) {6 ^ () 

This is enough to compute uf^ k with Eqn. 3-43, by recursion on 2g + k + n. We can 
also find a residue formula, based on the following principle. The divergent part of a 
rational function / at a pole A is given by: 

dz f(z) 

[f(zo)]-,x = Res — T -¥l (3-48) 

z^X Zq — Z 

u^ k (z ,I) is the sum of its divergent part near z = ±1, z = 1/zi, z = 1/sj and 
z = 0. After (vi), u^ ,k (z Q , I) has only a simple pole at z = 0, and it must be such that 
u^ k {z J) G 0(1/ z^) when z -> 0. Thus: 

To compute the residues at z — ±1, the term — d Zi (uj^ 1 (I)/ ■ ■ ■ ) do not contribute 
and we can replace S^ k by: 

£°< k (z,I) = {cu 9 n -\(z,z,I)+ 

+ E ^ x {z,J)u^- k '{z,I\J) 

^C/,0<9'<9,0</c'<fe 

+d z {^-\zJ))- X -^^ n k -\zJ) 

X \Z) 

(3-50) 

s ^>'> = ^) (3 " 51) 

3.6.3 Improvement for Co^'° 

This algorithm can be further simplified in the case k = 0, i.e the limit of a hermitian 
matrix model. Notice first that is obtained by the former algorithm without 
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computing any u^ ,k (z,I) with k 7^ 0. Then, one can prove the following "mirror 
relation" by a recursion using the loop equations: 



ouf(l/z,I) = z*u;f(z,I) 



(3-52) 



There is no similar mirror relation for u^' h (z,I) with k 7^ because of the term 
d z co^ k ~ 1 (z, I) in the loop equations. As a consequence, u^'° can only have pole at 
z = ±1. Indeed, no pole at z = Zi (resp. z = Sj) imply no pole at z = (resp. 
z = 1/sj) by the mirror relation. All the same, G 0(1/ z 2 ) when z — > 00 imply 
that u^'° is regular at z — 0. Besides: 



Res 



- Res dM'°(£,/) = 



(3-53) 



Therefore, we can compute u^'° in a simple way if the w^,'°'s with 2g' + n' < 2g + n 
are known: o;^' is the divergent part in the Laurent expansion of S^ ,k (z, I) at z — > ±1. 
We can again replace <S^'°(z, J) by <S^ ,0 (z, I) and write: 

dz 



u;f(z,I) 



Res 

z-»l,-l Zq 



(3-54) 



3.6.4 Standard form of the topological recursion 



Formula 3-54 is valid for one-cut solution of loop equations: we used the fact that 
u!^' is a rational function of its variables. It is false for other cases. However, we are 
going to transform slightly this formula into the form of the "residue formula of the 
topological recursion" [20], which still hold in more general settings. 

With the mirror relation, one can prove that S^'°(l/z, I) = z 2 (or use directly 
Eqn. 3-43 and Eqn. 3-44). Hence: 



uf(z,I) 



- Res dz ( — 

2 z-¥l,-l \Zq — Z Zq — I/Z 



Sf(z,I) 



dz 



Res 

z-n,-i 2y(z)x'(z) 



1 
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Zq - l/z 



9-1,0/ 



JCI,0<h<g 



(x(z) - x(zi)y 



Let us define: 



&2°{Z1, Z2) = ^2°{Zl,Z 2 ) + ~ 

n w Q /{z u z % Y = cof( Zl ,z 2 ) + 



1 dx(zi) dx(z 2 ) 



2{x{z l )-x{z 2 )) 2 

dx{z\) dx(z 2 ) 
{x{zi) - x{z 2 )) 2 



(3-55) 

(3-56) 
(3-57) 
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and if (n,g) ^ (2,0): 

wf( Zl , . . . , z n ) = "wf( Zl , . . . , z n y = u%°( Zl , ...,z n ) (3-58) 

We also define the point z ^ z such that x(z) = x(z). In general, it is defined at least 
locally around the simple zeroes of x' (the branchpoints at which we take the residues). 
Here, we have ~z= 1/z, which happens to be defined globally. 
Using the mirror property, we can rewrite: 

uf(z ,I) = ResdzK(z ,z)^w 9 n - + \'°(z,z,ir + 

Z— — 1 

JCI,0<h<g 

(3-59) 

This is the standard form of the topological recursion associated to a spectral curve 
(£, x, y), also called "residue formula". The residues are computed at the simple zeroes 
of dx on C. n za^°{z\, z%) n is the heat kernel of C, also called Bergman kernel. K(z, zq) 
is the recursion kernel. Here: 

K(z , z) 



Ay(z)x'(l/z) \z - z z -l/z 
1 ( \ 1 



2(y(^) - y(z))x'(z) \z -z Zq- z 
f*dzi»m 2 °(z>,z y 



(3-60) 



2{y{z) - y{z))x'{z) 
The last line is the expression of the recursion kernel in the general case. 

The residue formula (Eqn. 3-59) can be seen as a definition of some quantities 
(for 2 — 2g — n < 0) from the data of the spectral curve (£,x,y). These quantities 
satisfy some loop equations. When this formalism is applied to a spectral curve coming 
from a matrix model, it gives the topological expansion of the correlators of the matrix 
model. As a matter of fact, the residue formula allow to express all the solutions of 
loop equations (non necessarily with 1-cut). The quantities defined by Eqn. 3-59 
enjoy interesting properties, and can be computed with a diagrammatic technique (see 
the review [20]). 

Let us mention that another residue formula and a diagrammatic technique also 
exists in general for all ui^ k (the case of hyperelliptic spectral curves was treated in 
[11])- 

3.6.5 Example 

We apply the algorithm described in Section 3.6.2 to find the 1/N corrections to the 
1-point correlator. 
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(n,g,k) = (1,1,0). £i°{z) = a;^ ^, z) has a symmetry when z \-± l/z. So, killing 
the pole at z — Sj also kills the pole at z = 1/sj. Thus, u{'°(z) is the singular 
part at z = ±1 of S[ 1>0) (z). We find: 

1 1 



u{ {z) 



I67M (a) (z - l) 2 
1 1 



+ 



-M{b) + 4rfM'{b) 1 



4 7 M(6) (z + l) 4 4 7 M(6) (2 + l) 3 167M(6) 2 (2 + l) 2 

(n, g, k) — (1, 0, 2). This is the first case for which there is no short cut. 



0,2/ \ 



167M(a) (z - l) 2 4 7 M(6) (z + 



1 1 / ^ 1 \ 1 



5M' (6) 

7M(6) 1 16 ' (sj + l)(sj + 1) J ' 4M(6) 2 



E 



1 



+ 



— IF — 



E r « , a. 



(z + 1) 2 

jM'(b) 



( Sj + l) 2 s_y + 1 ; M(b) 



1 1 



Z+ 1 z 



E 



7 (M'(^)) 2 ( Sj + iy 



+ 



2s i 



1 



. F — 

jM'iXj) ( Sj - l)( Sj + 1)3 \ ( Sj - l)( Sj + 1) ^ s t - Sj 



(z-i/ Sj y 



1 1 



6 1 M'{x ] f 
M"( Xj ) 



(2M"'(x j )M'(x j )-3(M"(x j ))^ 



Sj + 1 



7 (M'(^)) 2 (s j + lf 
1 »2 



i-E 



si - Sj 



°j -^ + ^ + 2^-3 



7 2 M'(^) ( Sj - l)( S j + I) 3 I ^ ( s ' - s j) 2 s < - s i ( s o - W'i + !) 



E 



2 2(-2sj + 3) 

(siSj - l) 2 S;Sj - 1 



E 



, (sisj - l)(si>Sj - 1) 



z — I/Sj z 



3.6.6 The Bergman-tau function, F 10 

Finding F 10 requires the integration of ojI' (x) with respect to the variations of the 
potential, and we have: 

dz u)2°{z,z) 



uj\'°{zq) = Res 



z-ti,-i z — z 2y(z)x'(z) 
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(3-61) 



A closed formula (when there is no hard edges) was first given by Ambj0rn, Chekhov 
and Makeenko [3], and it was generalized to hyperelliptic curves by Akemann [2]. As we 
have seen, id^°(z, z) is closely related to the Bergman kernel of the spectral curve: the 
general method of integration of was studied in [21], and the solution was expressed 
in term of the "Bergman tau-function" . It depends on the positions of the edges (hard 
or not), and is defined as the result of integration of a closed 1-form. When there exist 
some hard edges, the result is slightly modified and was studied by Chekhov 4 [10]. We 
apply it in the one-cut case: 




(3-62) 



3.6.7 The Polyakov anomaly, F 2 

This term was first considered by Wiegmann and Zabrodin (Eqn. 5.20 of [28]). We 
adapt the expression given in [11] when there are some hard edges, in the one-cut 
case. By comparing with d a F 0,2 computed for the gaussian potential, we found that 
the result of [11] must be multiplied by 1/4 to be correct. Namely, the correct formula 
in our case is: 



^0,2 



k,k H In 

GO ^ 12 



M(a) 3 
M{b) 



a — b 
t 



(3-63) 



where: 



^0,2 



(I) i67r 2 



POi) 



Z2 — Z\ Z2 — 1 1 Z\ 



H P {z 2 )/t) 



C cxt is a contour outside the unit disk, such that indc oxt (z) = 1 whenever \z\ < 1, and 



ind 



■1 for all 1 < j < d- 1. 



Let us compute, for \z±\ > 1: 



A( Zl ) 



7T J C 



dz 2 



Z 2 -Zi Z 2 -1/Z! 



H P {z 2 )/t) 



7T JO I 9 

v (-i-E^ 008 ^)! 111 ^)/' 



k> 



^ kzf 



We have already computed this kind of integral in Eqn. 3-33. We just have to substitute: 

kv k — > 2j z\ 



Beware, the free energy is defined as — In Z in Chekhov, while we use the opposite convention. 
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in the result Eqn. 3-34. Therefore: 
A{z x ) = 

Then: 



d-l 



'd-l 



2lit 



8j\ I - 



1 fc>i 1 \i=i i 



k odd 



^0,2 



1 



1 



dp(*) 



/ **\—-—+t:(- 



The contour integral can be computed for each term of A(z). Individually, each term 
gives a meromorphic function to integrate. So, for each term, we can move the contour 
to surround the poles outside the unit disk, that is z = oo, and z = Sj (1 < Sj < d— 1). 
We obtain: 



^0,2 
(I) 



3=1 



3=1 



1=1 



d-l 



d-l 



^ — v ^ — v 2 1 ^ — r 2 1/ ^ — > 1 

+ 1^ k^k ~ l^k^k I l^^k 

3=1 fcodd 3 k>l 3 \l=l 1 



d-l 



d-l 

E 1 * 



In' S ' 



— N )+ 2 E ln f 1 -— ) 



Li=i 



s 1 



( Sj -i)( Sj + iy 



l<j<Kd-l 



We also write the second term of F 0,2 in terms of s,'s: 



12 V * 2 M W 



12 V t 



3=1 



i ( Sj - iy 



Therefore: 




(3-64) 

3.7 Unstable F^ k 

3.7.1 Computation of d a F 9 ' k 

We can repeat the argument which led to the residue representation for d a F 0,1 . It 
shows that the derivative term d z wf' k (z) has a vanishing contribution to the residue, 
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and we have 

Q aF9 ,k = 1 



Res -2y(z)uf(z) + £ 



a'.ife'/ \ g-g'k-k 1 



0<g'<9,0<k'<k 



x'(z) 



x'(z) 



(3-65) 

Another method is to extract directly cf' = —d a F 9,k from 7£f' , which is a by- 
product of the algorithm. Indeed, after Eqn. 3-42 and recalling y(z) ~ —r/(z — 1) 
when z — >■ 1: 



+ 0(z-l) 



(3-66) 



Then, after Eqns. 3-45-3-46: 



= _ r J- s _i+i Res de § f * (0 

-f-f Sj - 1 ^ 



(3-67) 



With this second method, we obtain: 



d a F^ 



d-1 

r x ^ 1 



-T 

7 2 M'(xj) ( Sj - l) 4 ( Sj - + l) 4 
d-l „2 

V r s j 

7 2 M '(^) ( Sj - + l)^ 

1 



* j -Sj + l 



(siSj-1) 2 SlSj - 1 (Sj - l)(Sj + 1) 



2s? - 3s? + 4s,- - 2 



+ 



, (SlSj - l)(si>Sj - 1) (Sj - l) 2 (Sj + 1) 



3.7.2 Algorithm for F 9 > k 

We state the so-called integration formula [11] for the stable free energies: 

F*> fe = — j Res dz<p(z)uf(z) (3-68) 

where is a primitive of y(z)x'(z): 

fZ 

ydx (3-69) 



For k — and g > 2, this formula simplifies into: 

F9 '° = 2~Yg z 5i G -i dZ ^K*^) ^ 3 - 7 °) 

This formula is not valid (or simply do not make sense) for unstable free energies. 
The correct formula were Eqns. 3-26-3-27 for F°'°, Eqn. 3-33 for F 0,1 , Eqn. 3-62 for 
F 1 ' , and Eqn. 3-64 for F°> 2 . 
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3.8 Normalization of the geometric quantities 
3.8.1 Position of the problem 

Let us define the partition function of a /3-matrix model by: 

Z N>p (a,V)= [ dA 1 ---dA A ,|A(A)| 2/3 exp(-^^y(A J ) | (3-71 

J]-ao,a]°° \ t i=l J 



It is known that, when the leading order (N — > oo) of W\(x) has only one cut in the 
complex plane, there exists an asymptotic expansion (N — > oo) of the form: 



\nZ Nt0 (a,P,V) = Dp\nN+ £ (^J^) (^"^) + 

(3-72) 

Dp depends only on 0, and d 9,k is a constant. They do not depend on a, V and N. 

To this matrix model, we associate a plane curve S a y, which do not depend on N 
and (3: 

S a y : y = <'°(x) - (3-73) 

Out of the curve were construct some geometric quantities F 9 ' k (S a y). We can form a 
formal asymptotic series in N: 

?NAS a y)=^(^) 9 (v^-^y ^\S ay ) (3-74) 

The theorem completely proved in [11] (and only sketched here) is that there exists 
constants d 9,k such that 

V<?, k > f 9 > k (a, V) = F 9 >\S a y) + d°> k (3-75) 

In other words, there exists a constant Cn,/3 depending only on N and 0, such that: 

Z Nt/3 (a, V) = C N £ exp (T N ^(S a y)) (3-76) 

and lnC^ has a topological expansion, modulo an extra ln(iV) term. 

If one wishes to test Eqn. 3-76 with Monte-Carlo computations of Z, or to compare 
the predictions for the statistics of the maximal eigenvalue to Tracy- Widom laws, it is 
important to know the normalization constant CV^. 

3.8.2 Value of C N ,p 

It can be computed explicitly by specializing to some a and V such that: 
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• one has a close formula for F 9,k for any g, k. 



• one has independently a close formula for the eigenvalue integral. 

Comparing the two gives the constant p. Such expressions are available for the 
Gaussian potential (V(x) = x 2 /2) in the limit a — > — oo. In this case indeed, the 
eigenvalue integral can be reduced to a Selberg integral on one hand, and unstable F 9)k 
vanish on the other hand. We do the computations in Appendix A, and the result is: 



Cjv,/3 = exp 



TV 2 In N + f3N 2 ( - - In /3 



{fi - l)JV(-l + ln(27r) +ln/3)] • 2?"^ 



nf = ir(i+j/3)) 2 
r(i + /?)^r(i + /3A^) 

(3-77) 



We shall see from this formula how to find the constant prefactor in the asymptotic 
of the (3 analog of Tracy- Widom laws for arbitrary (3. 

4 Example: Gaussian /3 ensemble 

We consider the Gaussian model, i.e with a quadratic potential V(x) = \: 

X Nj p(a)= f dA 1 ---dA JV |A(A)| 2/3 e-^ E -i A ? (4-1) 

J]-oo,a] N 

We also write X N ^ for the unbounded eigenvalue integral X N ^(oo). Let us note F 9 ' k (a) 
the geometric quantities of this model, and write: 

F(a) = J2 v 2 ' 29 hk F9 ' k (a) (4-2) 

9,k 

We are interested in: 

V(a) = P[A max < a] = = ?M (4-3) 

Cn,/3 = e Djv ' 8 is the normalization constant relating the series constructed out of the 
geometric quantities to the genuine eigenvalue integral. Its value is given in Eqn. 3-77. 

Notice that Xpf,p(a) is a function of a, t, but we have a priori a relation of homo- 
geneity: 

X N ,p(a) = \ar 2+ ^ N I(t/\af) (4-4) 



where: 



J(r)= / dA 1 ---dA 7V |A(A)| 2/3 e-^ E -i A ? (4-5) 

J]-oo.±l] N 



-oo,±l 
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Thus: 



F°'°(a) = t 2 In \a\ + function of t/\a\ 2 (4-6) 

F 0,1 (a) = -tin \a\ + function of t/\a\ 2 (4-7) 

elseF 9 ' fc = function of t/ 1 a| 2 (4-8) 

4.1 Theoretical results for arbitrary a 
4.1.1 Density of eigenvalues 

The eigenvalue support is a single interval 5 if a < 2yfi. In this case, it is of the form 

[b, a] where b = | (a — 2\Ja? + 12t) . The parameters a = 9 ^ and 7 = ^ are thus 
given by: 

a = \{2a - y/ a? + I2t) (4-9) 



7 = -(a + VaF+m) (4-10) 
o 



and we have: 



a 2 n 7 2 

w = — + 7 u 1 = "7 ^2 = y (4-11) 

r = —a u = a u\ — 7 (4-12) 
The density of eigenvalue is: 



1 x — b 

p(x) = M(x)\ for 16 [6, a] 

2ix \ a — x 

M(x) = x - 27 

When a = 2y/t, the eigenvalue support is [b, a] = [—2y/i, 2y/t] and we recover the 
Wigner semi-circle law (with t — 1/2 and a = \/2 usually) : 



= — V 4t — x 2 
2n 

For a < 2\/i, p(x) has another zero on the real line at x — X\ — 27, which corresponds 

to S — S\\ 

si = — (2-1 - a + y/a{a - 47) j 

defined by |si| > 1 and x(si) = x±, where x(z) = a + 7 (z + |J is the Zhukovsky map. 
For a < we have s G [1, oof. 



5 We will not analyze the case where a > 2\/£, which cannot be reached naively as below 
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4.1.2 First terms in F(a) 

The formulas for the nondecaying terms in the free energy yield: 

F°.o W = «¥ + t- t (y* + ?P) + i']n'y (4-13) 

^» = -S)-(- + t) + (t^ 2 )j 

F 1 ' ^) = -^ln[l6r 4 7 4 a 3 («-47)] (4-15) 

F o,2 (a) = 1 ln f 16t -2 7 4 g^-Zi)! ^ (4-16) 



12 V + 



Indeed, the dominant term (Eqn. 4-13) was obtained in [14, 15] for all /3 by Coulomb 
gas techniques. 

Their derivatives with respect to a are maybe easier to handle: 

d a F°'°{a) = 7 a 2 (4-17) 

d a F°'\a) = + -4 7 )) (4-18) 

^(a) = - ^ (4-19) 
a(a — 47)^ 

For example, in terms of a and t only: 

F^(2Vt)-F '°(a) = ^a 2 -^a 4 --^-(a 2 + 30t)v^+12i 

3 zlo zlo 



-t 2 ln 



+ Va 2 + 12t 
6^ 



(4-21) 



iri/ \ 7 ln 3 ln 2 ln t 1 , . 9 

F ' (a) = + lna 2 + 12t 

v 7 24 8 6 48 V 7 

-—In [a (-a 6 + 6a 4 t + 36a¥ - 432t 3 ) + 
24 ' 



(-a 6 + 12a 4 t - 54a 2 t 2 + 216t 3 ) Va 2 + 12t 



(4-22) 
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4.1.3 Example: a = 

V(0) = y^£ F ^ is the probability that all the eigenvalues are negative. The first 



terms in F(0) are 



F°'°(0) 
^0,1(0) 

F x '°(0) 
F°' 2 (0) 



3 ln3 lnt 



S + -( 1 + £)-N«-T 

In 2 In 3 



31n3 1, / 2 
h - In 1 + — = 

8 2 V \/3 



(4-23) 

(4-24) 

(4-25) 
(4-26) 



4.2 Numerical simulations: principles 



To verify the analytical predictions derived in the previous sections, we simulated the 
joint distribution of eigenvalues in Eqn. 1-2 in the case of the Gaussian ensemble: 



P(X 1 ,...,X N ) = B M>N l[e-^^ ] J |A, - A 



|2rf 



i<j 



— Bm,n e 

where the effective energy E [Aj] is given by 



N \ 2 

i=l 



(4-27) 



(4-28) 



i<j 



In the simulations, we fix t = 1/2. We sample this probability distribution using a 
Monte Carlo Metropolis algorithm (see [29]). 

4.2.1 Standard Metropolis algorithm 

We start with an initial configuration of the Aj's taken at random within the range 
[-2y/i,2y/i\ = [-V2,V2] (the average density support). At each step, a small mod- 
ification {Aj} — > {A'j} is proposed in the configuration space. In our algorithm, the 
proposed move consists of picking at random an eigenvalue Xj and proposing to modify 
it as Xj — > Xj + e where e is a real number drawn from a Gaussian distribution with 
mean zero and with a variance that is set to achieve an average rejection rate 1/2 (the 
proposed move is indeed accepted or rejected as explained below). 
The move is accepted with the standard probability 

P(X 1 , . . . , X' N ) 



v 



mm 



P(A 1 ,...,A 



mm e 



(4-29) 
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and rejected with probability 1— p. This dynamics enforces detailed balance and ensures 
that at long times the algorithm reaches thermal equilibrium (at inverse "temperature" 
(3) with the correct Boltzmann weight e~@ E ^ Xi ^. 

At long times (from about 10 6 steps in our case), the Metropolis algorithm thus 
generates samples of {Aj} drawn from the joint distribution in Eqn. 4-27. We can then 
start to construct a histogram for the maximal eigenvalue (by keeping the value of A max 
about every 100 steps). This histogram describes the probability density function of 
Amax, which can be compared to our analytical predictions: P(A max = a) = d a V(a) with 
V(a) given in Eqn. C.13 for large N. More precisely, we plot here the rate function 
$(a) = — lnP(A max = a). We compare the numerical data $ num (a) with analytical 
results $(a). 

However, as the distribution of the maximal eigenvalue is highly peaked around its 
average, a standard Metropolis algorithm does not allow to explore in a "reasonable" 
time a wide range of values of A max . The probability to reach a value A max = a decreases 
rapidly with iV as e - N2 ®2( a )+°( N ) w h ere $ 2 (a) is positive (for a different from the mean 
value). Therefore, we modified the algorithm in order to explore the full distribution 
of A max and to compare it with our analytical predictions for the (left) large deviation 
of the distribution of A max . 

4.2.2 Modified algorithm: conditional probabilities 

It is difficult to reach small (or large) values of A max with a standard Metropolis algo- 
rithm. The idea is thus to force the algorithm to explore a region A max < a c for different 
values of a c . We thus add in the algorithm the constraint A max < a c . More precisely, 
we start with an initial configuration that satisfies the constraint A max < a c . At each 
step, the proposed move {Aj} — > {A^} is rejected if A' max > a c . If A max < a c , then the 
move is accepted or rejected exactly with the same Metropolis rules as before. Because 
of the new constraint A max < a c , the moves are rejected much more often than before. 
Therefore the variance of the Gaussian distribution P(e) has to be taken smaller to 
achieve a rejection rate 1/2. 

We run the program for several values of a c (about 60 different values on the left of 
the mean value) and we construct a histogram of A max for each value a c . This gives the 
conditional probability distribution P (A max = a|A max < a c ). Again, as the distribution 
of the maximal eigenvalue is highly peaked, the algorithm can only explore a very small 
range of values of a - even for a large running time (about 10 9 steps). The difference 
with the previous algorithm is that we can now explore small regions of the form 
a c — r] < a < a c for every a c , whereas before we could only explore the neighborhood 
of the mean value. 
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The distribution of the maximal eigenvalue is given by 

P (A max = a) = P (A max = a\ A max < a c ) ■ P (A max < a c ) (for a < a c ) . (4-30) 

Therefore the rate function reads: 

- lnP(A max = a) = - [lnP (A max = a|A max < a c ) + InP (A max < a c )] . (4-31) 

The histogram constructed by the algorithm with the constraint A max < a c is the rate 
function "^(a) = — InP (A max = a|A max < a c ). It differs from the exact rate function 
$(a) = — lnP(A max = a) by an additive constant that depends on a c . In order to 
get rid of this constant, we construct from the histogram giving $ ac ( a ) the derivative 
d# a e (°) _ T^jg derivative is equal to d ^ Q - ) and the constants disappear. 

4.3 Numerical simulations: results 

We now come back to $(a) from its derivative using an interpolation of the data 
for the derivative and a numerical integration of the interpolation. This allows to 
compare directly the numerical results $ num (a) with the theoretical rate function 
$(a) = — lnP(A max = a) = —\n[d a V(a)]. The interpolation and numerical integra- 
tion also help to get rid of the fluctuations in numerical data (by smoothening the 
plot) and thus increase the precision of the results. 

We can follow the same steps to explore the region on the right of the mean value 
by adding in the simulations the condition A max > a c (instead of A max < a c ) for several 
values of a c . At the end, we can reconstruct the full plot of the probability density 
function (pdf) of A max . There remains only an additive constant in the rate function 
that is not yet determined (as we perform a numerical integration of the derivative). 
This constant can be numerically determined by imposing that the integral (numerical 
integration of the data) of the pdf of A max is one (normalization of the probability). 
We can now compare the simulations with analytical predictions. 

Our analytical predictions give (see Eqn. C.13): 



V(a) 




exp {F(a) + lnC N> p - lnX^} 



(4-32) 



with 




ln(2vr) In (3 



+ Kp 



(4-33) 



2 2 
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and F(a) is given up to order one (for large N) by 

TV 2 /V 

F(a) = —f3F 0fi (a) + —IB - l)F°Ha) + F lfi (a) + (8 + B~ l - 2)F°' 2 (a) (4-34) 
t 2 t 



Thus 



where 



$(a) = -lnP(A max = a) = -\n[d a P(a)) 

= - {F(a) + hxC N>p - \nl N> p + In [d a F(a)}} 

= N 2 $ 2 (a) + N<S> 1 (a)+\nN(j) + %(a) + --- (4-35) 



* 3 (a) = -I/H*°(a)-/?(2-^ 

*i(<0 = -^(/3-l)F°' 1 (a)-(/3-l)(-^ + ln(2 7 r) + ^ 

3 + 8^-3 
12 

$ (a) = -F 1 ' (a)-(/3 + /3- 1 -2)F°' 2 (a)-ln(9 a F '°(a)) 

5-2(^)X hn>) M _^_ ln f4 
6 y 2 2 w 2 

We typically run the simulations for N = 30 with 10 9 iterations. 

Fig. 1 shows a plot of the rate function $ num (a)/iV 2 = — lnP(A max = a)/N 2 for 
different values of 8: in red, /? = | = 0.625 ; in blue, 8 = 1 (Gaussian Unitary 
Ensemble) ; and in orange, 8 = | = 1.6. Numerical data $„„(o) (disks, squares and 
triangles) are compared with analytical predictions for the left tail $(a) = A^ 2 $2(o) + 
A r $i(a) + (ln A^)0+$o(o) (solid lines). In this figure, we see that the agreement between 
numerics and theory is very good, at least to dominant order for large N (order iV 2 ). 

Fig. 2 shows the next order in the large iV expansion (order N) and figure 3 shows 
the order one for large N. Up to order O(N), the agreement between numerics and 
theory is very good (see Fig. 1 and Fig. 2). To order 0(1), the theoretical lines do not fit 
numerical data as well as to order O(N) or 0(N 2 ) but the agreement is still quite good. 
At this order, the number of significant digits that is required is very high. A better 
precision would require a longer running-time (to have more samples) and a refined 
histogram. Also, the determination of the constant depends on the reconstruction of 
the full histogram (both left and right tail). However, we can not go numerically toward 
infinity, we have to choose a finite interval. The precision could thus be increased by 
exploring a wider range of values (by going deeper in the left and right tails). 
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Figure 1: Case of a quadratic potential (Gaussian ensemble): rate function Q(a)/N 2 = 
— lnP(A max = a) /N 2 = —\n[d a V(a)} /N 2 plotted against a for different values of (3: 
(3 = 5/8 = 0.625 (red disks and red line), (3 = 1 (blue squares and blue line) and 
/? = 8/5 = 1.6 (orange triangles and orange line). The points (disks, squares, triangles) 
are numerical data obtained with the modified Metropolis algorithm. The solid lines 
are the analytical predictions for the left tail of the large deviation of the pdf of A max : 
V(a) is given in Eqn. C.13. 
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Figure 2: Case of a quadratic potential (Gaussian ensemble): order N of the rate 
function, ($ num (a) — A r2 $ 2 (a))/A r , plotted against a for different values of (3: (3 = 
5/8 = 0.625 (red), (3 = 1 (blue) and (3 = 8/5 = 1.6 (orange). The points (disks, squares, 
triangles) are numerical data obtained with the modified Metropolis algorithm. The 
solid lines are the analytical predictions: (A r $ 1 (a) + (In A?") <fi + <&o(a))/N. 
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$(a) - iV 2 $ 2 (a) - N$i(a) - In N <fi 




Figure 3: Case of a quadratic potential (Gaussian ensemble): order one of the rate 
function, $ num (a) — N 2 & 2 (a) — AT$ 1 (a) — (In N) 0, plotted against a for different values 
of (3: (3 = 5/8 = 0.625 (red), (3 = 1 (blue) and = 8/5 = 1.6 (orange). The points 
(disks, squares, triangles) are numerical data obtained with the modified Metropolis 
algorithm. The solid lines are the analytical predictions: $o( a )- 
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5 Scaling limit at the edge of the spectrum 



Now, we wish to take the scaling limit N — > oo, a — y a the edge of the spectrum, in 
a one-cut model with arbitrary polynomial potential. Let us mention that using an 
asymptotic expansion of eigenvalue statistics to study limiting distribution is not a new 
idea. For instance, within GUE, Gustavsson [22] has found limiting distributions for 
the fc-th eigenvalue in various regimes, using an asymptotic expansion for the first few 
terms of eigenvalue correlation densities. 



5.1 Blow up of the spectral curve and unstable correlators 



0,0 



0,1 



Formula 3-68 computes the stable F 9 ' k (a) only in terms of the spectral curve y, co 2 
and UJi 

( a) ' 1 viz) = m n*zl ( } 

(5 " 2) 



U2°(zi, z 2 ) = — : — (5-3) 

(ziz 2 - V 2 



We shall use the notation F 9,k (a) = J-* flj fe(E a ) for the geometric quantities associated 
to the curve E a . 

Let us define e such that si — 1 ~ y/e, and e — > when a — > a*. Let us also assume 
that 7 — > 7 7^ and Sj — > s) ^ 1 for 2 < j < d — 1. With the rescaling z = 1 + ^/ecr, 
we have: 

, v x / x(z) = (a - 2e7) + 7e x(cr) 

~ 1 »W = A7 + ^ y(a) < 5 " 4 > 
^•'(z) ~ e~ 1 ' 2 S°'V) (5-5) 

0,0/ \ -1 -~0,0/ \ ir a \ 

lo 2 ' (zi, z 2 ) ~ e u 2 '(cri,cr2j (5-6) 

d— 1 

i=2 



with the constant: 



We have: 
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Legend • branchpoints (zeroes of cLr) 

unphysical sheet 
Sj or 1/sj (zeroes of y) 



Figure 4: When a moving singularity {x = x(sj) here) merges to a fixed singularity 
(the hard edge at x = a), the blow up around a of the spectral curve (here, of the 
form Eqn. 5-1) can be generically described by the situation on the left. Indeed, 
(£) = (C, x, y) is a universal object. We may call it the Tracy- Widom curve. In our 
parametrization, the physical sheet correspond to Re a > 0. There is only one simple 
branchpoint, at o = 0, and we have a = —o (which is globally defined), y has a pair 
of zeroes at a = —1,1. 



In other words, E is the plane curve of equation y 2 = x + 1/x — 2. 
5.2 Scaling limit of unstable F 9,k 

It is easy to obtain the scaling limit of unstable free energies from their expressions. 
We find: 

F°'°(a) = F°' (a) + ((7M + ) 2 / 3 e) 3 F '° + o(e 3 ) (5-10) 
F o,i (a) = F°' 1 (a) + ((7M + ) 2 / 3 e) 3/2 F°' 1 + o(e 3 / 2 ) (5-11) 
F lfi (a) = ln((7M + ) 2 / 3 e)F 1 '° + 5 1 '° + o(l) (5-12) 
F°< 2 (a) = ln((7M + ) 2 / 3 e)F°< 2 + <5 ' 2 + o(l) (5-13) 



12 



with: 



and: 



= -|, (5-14) 
F ' 1 = ^ (5-15) 
F 1,0 = -I, F°> 2 = - (5-16) 



5 1 ' = -— ln(2 10 M+M_7V 4 ) (5-17) 



«5°< 2 = -In 
12 



24 



2<j<l<d-l 



F°'°, F 0,1 , F 1 ' and F ' 2 are universal and can be considered as the unstable free 
energies associated to the curve S. On the contrary, Sip and 5 0,2 are non universal. 
They are relevant to compute the constant prefactor in the scaling limit of F(a). 

5.3 Scaling limit of stable F 9,k 

The scaling of the correlators is very simple to obtain from the residue formula 
(Eqn. 3-49). The scaling of unstable correlators are computed at hand, and the re- 
sult for stable correlators is proved by recursion, cu®' is special: 

uj°Az) = \v\a) e - 7M + ^ W) + 0(e 2 ) 
And, if (n t g,k)^ (1,0,0): 

U% k (Zi, ...,Z n ) = (7M + e 3/2 } 2-2 9 - fc -n Q g, k ^ ^ ^ + ( g_ lg) 

The Q^ h are the result of the residue formula applied to £ and the data of to®' and cD'j'' 1 . 
The standard form of the topological recursion (Formula 3-59) computes with the 
limit of the recursion kernel: 

*o) = \ - f 2 t4 2 r (5-20) 

4cr(cH - l)(cr 2 - cr ) 

In the stable the leading behavior correspond to diverging power of e. Thus, only 
the poles at infinitesimal distance of the hard edge are relevant to find 0^ h . In other 
words, the poles at Zi — 1 are replaced by poles at a,i = 0, the poles at Zi — 1/si are 
replaced by poles at <7j = —1, the poles at Zi = are replaced by poles at Oi = —a h 
whereas the poles at z = and z — 1/sj do not give contributions to u)^ k . 
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The integration formula for stable F 9,k (Eqn. 3-68) yields the scaling: 

F 9 ' k (a) = (jM + e 3/2 ) 2 ' 29 ~ k F 9 ' k 
where F 9,k are the geometric quantities associated to the curve E: 



J?9r 



1 



2 - 2g - k e-w-i 
is the following primitive of y dx: 

[a) = I ydx 



Res da 4>(<j)u)f' k (a 



(5-21) 
(5-22) 

(5-23) 
(5-24) 



5.4 Summary: double scaling limit 

To find the scaling limit of F(a) when a goes to the edge of the spectrum, a scaling 
variable appears naturally: 




(5-25) 

Notice that the N 2 ^ 3 scaling is the same for all /3, as pointed out long ago, and proved 
recently in [24]. When t] is of order 1, each term in the topological expansion of F(a(r))) 
is also of order 1. It is the so-called "double scaling limit": 

2-2g-k 

- I /3 1 - s (l-/3- 1 ) fc F^(S a ) 



£ (t 

g,k>0 v 



In 



+6 1 ' + (P + /T 1 - 2)5 ' 2 + ^ ^(l-r 1 )^ 1 ^ 29 ^^ (5-26) 



g,k>0 



where it is implied that 77° should be replaced by \n(r]). In particular, for hermitian 
matrices (/3 = 1): 



F(a) 



t 2 12 V * 



+Si,o + J2 r i 3[1 ~ g) F9,c 

9>0 



(5-27) 



44 



The curve £ was described in Section 5.1. We call it the curve the Tracy- Widom curve, 
because its geometric quantities F 9,k are related to the coefficients of the asymptotic 
expansion of the left tails of the Tracy- Widom laws (see Section 5.6.2). 

Eventually, the complete answer for V(a) = P(A max < a) is: 

oo 

1-fl (l - /3 _1 ) fc rfi {2 - 29 ~ k) F 9 ' k (5-28) 

g,k>0 

Notice a posteriori that a must be such that the divergent terms when N — > oo cancel 
in this expression. 

5.5 First few F^ k 

We give a few of these geometric quantities associated to S. We have already seen: 

- 9 
F°'° = -- 

3 

- 4 
F - 1 = - 

3 

- 1-1 



12' 12 



We can also compute: 



£1,1 _ 31 p 0:3 _ _9_ 

2 7 -3' 2 6 
£2,0 __ £1,2 _ 487 90,4 _ 595 



2 9 ' 2 12 -3' 2 n -3 

~ 2>1 _8831_ ~ 1;3 9281 ~ 0)5 19977 

2 i7. 3 ' 2 17 -3 2 ' 2 16 -3 

£3,0 = 63 ^ 2 , 2 = 62761 p M 539193 p , 6 577879 
2 14 ' 2 20 ' 2 23 ' 2 22 

5.6 Example : Gaussian j3 ensemble 

To illustrate the general framework explained above, we carry on direct computation 
in the Gaussian case from the results of Section 4. The edge of the spectrum is located 
at a = 2\ft. This value is actually the mean value of the largest eigenvalue of the 
Gaussian model: 

lim E(A max ) = 2Vi (5-29) 

N— >oo 
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F(a) gives the large deviation of the distribution of A max on the left of the mean 
value, which describes atypical large fluctuations of order 1. On the other hand, it is 
known that the typical small fluctuations of order N~ 2 ^ 3 around the mean value are 
described by the Tracy- Widom distributions. We expect that the two regimes match 
smoothly and thus that a scaling limit of F(a) when a — > (2y/t) coincides with the 
limit £ — > — oo of the Tracy- Widom distribution F2 J g(£)- More precisely, we write: 

a = 2y/i(l - e) (5-30) 

and the correspondence 

e = _ C/3 TV" 2 / 3 £ (5_31) 
is expected for some constant > 0. 

5.6.1 First terms in F(a) 

For bookkeeping, let us write several useful expansion when e — > 0: 



a 



1 



e e 2 e 3 



VMl + - + - + -) + (e S) (5-32) 



e e 2 e 3 e 4 



Vt[l-- + — + — + — - + 0(e 5 ) (5-33) 
2 16 64 1024 ' 



The other zero of p(z) outside the unit disk is: 

r e 7e 3 / 2 5e 2 103e 5 / 2 9e 3 
31 = 1 + ^ + 2 + ^6- + 16 64 

615e 7 / 2 45e 4 10187e 9 / 2 5 . . OA . 

+ W + K)24 + ^2W + 0(e) (5 " 34) 

We see that the definition of e through Eqn. 5-30 is compatible with the notation in 

Section 5.1. Moreover, the natural scaling variable is r] = t~ l / 2 e. We have for the 
nondecaying terms of F(a): 

F^(2Vt) = t\- 3 -+ l f) (5-35) 

F°>\2Vi) = *Q-ln(27r)-^ (5-36) 

F 1,0 (2 y/i) diverges logarithmicallly 
F 0,2 (2y/t) diverges logarithmically 
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And: 



F°>°(a) - F°>°(2Vi) = t 2 (-^- + ^j+0(e 5 ) (5-37) 

/A J/2 2 5/2 3 

F°' 1 (a)- J P°' 1 (2Vt) = M- 

y ' y 1 \ 3 2 20 24 

23e 7 / 2 e 4 i43e 9 / 2 \ 

+ " 7^ + "77777 + 0(e 5 ) (5-38) 



(5-39) 



F o, 2( 



896 128 18432 / 
n 2 5e 5e 2 
T~ + 64 + 1024 

lne 71n2 13e 
~8 12~ + ~2 ~ 64" 

13e 3 / 2 69e 2 147e 5 / 2 97e 3 
+ _ 96 1024 + 5120 ~ 6144 

103e 7 / 2 83e 4 35005e 9 / 2 ^. 5 . 

1 h 0(e 5 ) (5-40) 

114688 131072 9437184 y J K 1 

5.6.2 Scaling regime of F(a) 

We have argued in Section 5 that, when e is of order N~ 2 ^ 3 , all the terms in F(a) 
become of order 1. This can be checked explicitly on Eqns. 5-37-5-40. N~ 2 / 3 is indeed 
the typical scale expected of small fluctuations around E(A max ). Let us write V(a) in 
the scaling variable r\ such that 

a = y/i(2 - 2N~ 2/3 r)) (5-41) 

and taking into account the first terms F 0,0 , F 1,0 , F ' 1 and F ' 2 . When we use the 
value in the large N limit of C^,p (Appendix C.3) and In,/3 (Appendix C.2), we find 
as expected that the following limit exists: 



lim V[a = 2Vi(l - N- 2 'S)} = V* G ^( V ) (5-42) 
For a better comparison to the literature, we set £ = — 2rj, and we find: 



^g/3e(^) = T H ex P 



(5-43) 

T/3 is a constant given by: 

/17 25, n „ u \ ln(2vr) In/3 

lnT/3= IT " 24 (/3 + /? O 1112 - ~2 ~^ + ^ (5 " 44) 
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Kp is the finite part in the large N asymptotic expansion of ^ J=1 lnr(l It is 

studied in Appendix C.l, and computed in a somewhat explicit way for rational values 
of j3 in Appendix C.4.2. 



Value of 



Hrp) 



Numerical value of 



+ C'(-i) 



In 2 
' 12 

371n2 , C (-1) 



48 ~ 


2 


97 In 2 


71n3 


72 


36 


87 In 2 


La(2ir) 


32 


4 ' 


11 In 2 | 


C'(-i) 


48 ~ t ~ 


2 


971n2 | 


51n3 


72 ' 


18 


185 In 2 


ln(27r) 


96 


4 



InjgTr) , lnr(l/3) , C'(-l) 
6 ' 3 ' 3 
lnr(l/4) , c'(-i) 

2 ~ r 4 



0.8723714 
0.5395545 
0.3071491 
0.1752911 



1/2 
1/3 
1/4 



ln(27r) 
+ 



lnr(l/3) , C'(-i) 
6 1 3 ' 3 
lnr(l/4) , c'(-i) 

2 < 4 



0.7854042 
0.5160081 
0.3034417 



2/3 
3/2 



65 In 2 _,_ 2 In 3 _|_ ln(27r) 



144 1 9 
17 In 2 _ 25 In 3 
72 



inr(i/3) , C ; (-i) 

Q ' 



144 



+ 



Infgjr) 



lnT(l/3) 
3 



(-D 



0.8882751 
0.7051367 



5.6.3 Comparison to Tracy- Widom law 

This expansion should match the left tail of Tracy- Widom law [25, 26] for the usual 
values of = 1/2, 1,2. In fact, the leading term in e 3 was obtained in [14, 15] from 
Eqn. 4-13. The correspondence between the scaling variable £ in the large deviation 
and the Tracy- Widom variable x should be the following: 

• Gaussian Unitary Ensemble (/3 = 1). 

F 2 (X) = P*gve(-x) (5-45) 
If we specialize Eqn. 5-43, we obtain: 

P 2 (X) = 2 1 /2V'C- 1 )^(l + 0(| X r 3 )) (5-46) 

X->~oo \X\ ' 

Notice that only F 9 ' are involved in Eqn. 5-28 when (3 = 1, so the 0(|x| _3//2 ) is 
in fact a 0(|x| -3 ) according to Eqn. 5-21. 

• Gaussian Orthogonal Ensemble (/3 = 1/2). 
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We obtain: 

Ixl 3 Ixi 3 / 2 

Pi(x) = 2- 11 /4B e a-D/2 g " 3 ^ (i + 0(| X |- 3/2 )) (5-48) 

• Gaussian Symplectic Ensemble (0 = 2). 

F 4 (X) = ^se(-2- 2/3 x) (5-49) 

We obtain: 



F 4 ( X ) = 2 -3 5 /48 e a-D/2 ^ 3 2 (l + 0(|x|- 3/2 )) (5-50) 

X— >— oo X 

These expansions match earlier results. Usually, the constant term is considered 
as the difficult part of the asymptotics. For = 1, this constant was first obtained 
by P. Deift, A. Its and I. Krasovsky by Riemann-Hilbert asymptotic analysis [16]. For 
— 1/2 and 2, it was determined by J. Baik, R. Buckingham and J. DiFranco by a 
technique allowing representation of total integrals of the Hastings-MacLeod solution 
of Painleve II [5]. In these works as in ours, an important step was the evaluation 
of the asymptotic of a Selberg-type integral (see Section 3.8 and Appendix A). We 
remark that for ^ 1/2, 1, 2, the choice of the normalization of the scaling variable is 
somewhat arbitrary, since at present, there is no other definition of a " Tracy- Widom" 
distribution than: 

r 2/3 - /^G/3E 

Our T/3 gives the constant of F 2l a(x) = Vq^t] = — x/2). 

6 Conclusion 

• We showed how to determine to all orders the left large deviation function of the 
maximal eigenvalue of a random matrix in the ensemble, for all values of 0. 
We have presented in detail the case where the large N spectrum is connected 
(one-cut case), but the method combined with existing literature would give as 
well the multi-cut case. We performed numerics to check our results on the case 
of the Gaussian model at orders 0(N 2 ) and O(N) with a very good agreement, 
and at order 0(1) with a less convincing agreement. 

• We described the repartition function F 2( g(x) of the maximal eigenvalue when 
it reaches the edge of the spectrum from the left (existence of this distribution 
was proved in [24] with probabilistic methods). More precisely, we gave its full 
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asymptotic expansion when x — )• — oo in terms of geometric quantities associated 
to a spectral curve of equation: 

£ : y 2 = x + --2 (6-1) 

x 

Since it is well known that F1/2, F 1; F 2 are the Tracy- Widom laws, we have 
related these Tracy- Widom laws defined from the world of Painleve equations and 
integrable systems, to the topological recursion defined from algebraic geometry. 
For /3 — 1, in a related work [8], two of us proved directly that the resummation 
of the symplectic invariants of S* are solution to an ODE related to Painleve 
II, as it should be according to the work of Tracy and Widom. For general (3, 
the topological recursion allows to define a F 2( g for any value of (3. Here, we 
have obtained the three first terms of its left tail, the constant term, and we 
may compute recursively the next terms. We notice that the constant term was 
computed in three parts, similarly to [5] where the problem was adressed for 
= 1/2,1,2. Eventually, it would be interesting to relate the F 9,k [S*], to the 
parabolic differential equation obtained recently in [6] from which the Tracy- 
Widom (3 distribution can be extracted. Nevertheless, at present, nothing is 
known about the relation between integrable systems and (3 matrix models. 

• From the mathematical point of view, our work is based on two assumptions : 
existence of a 1/N expansion, and possibility to take a double scaling limit. The 
first one will be justified in a forecoming work [9], whereas the second one seems 
hard to justify rigorously at present in absence of integrability. If one accepts 
the 1/N expansion, the loop equations imposes the special /3-dependance of the 
expansion. From this, one observe the following duality (3 ■H- 1/(3 at the level of 
left-tail asymptotics: 

^e(x) = /^7W/3 2/3 x) (6-2) 

where ~ is obtained by taking the other branch of the square root (all half- 
integer powers of s of the expansion appear with a minus sign). Such a duality 
is readily observed on the known results for F± and F 4 (compare Eqn. 5-48 and 
Eqn. 5-50). Whether this duality has a counterpart for the functions themselves 
is an open question. 

• One could ask for a combinatorial interpretation of the symplectic invariants 
of the Tracy- Widom curve J 73 ' ^) (let alone T 9 ' k (S)) for fixed g. A related 
challenge would be to have a closed formula, at g fixed, for F 9,0 and u 9 ; of the 
Tracy- Widom curve. 
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Appendix 



A Determination of the normalization C^,p 

When a — > — oo, we have the asymptotics: 

F°'°(a) = -^-t 2 ln|a|-^ + t 2 lnt + 0(|a|- 2 ) 
F°'\a) = t\n\a\ + t - tint - t\n(2n) + 0(\a\~ 2 ) 
F lfi (a) = -^ + 0(|a|- 4 ) 

F°< 2 (a) = ^ + 0(|a|- 2 ) 

and one can prove that all other F 9 ' k are 0(\a\~ 2 ). As F(a) = J2 g k>0 v 2 ~ 2g ?i k F g,k (a), 
we have the following result: 

lim \F(a) - ((3N 2 {-\a\ 2 /(2t) -lnlal) + (B - l)Aln \a\)] 
= f3N 2 ( ~ + In t] + (/3 - 1)N (1 - ln(27r) - In t) 



2 

+ 5 
3 6 



In 2 



(A.l) 



On the other hand, consider the bounded eigenvalue integral 

X NtP {a)= [ dA 1 ---dA JV |A(A)| 2/3 e-^ E -i A ? (A.2) 

J]-oo,a] N 

One knows that: 

T N ,p{a) = C NiP e F ^ (A.3) 

We want to find Cn,p by comparing the result of Eqn. A.l to the large N asymptotic 
of X/v i/ g(a). Heuristically, when a — > — 00, we do not see a Gaussian weight, but the tail 
of it starting from the point a, which looks like a decreasing exponential. Precisely, for 
any finite N, with the change of variable: 

Mi = ^M (a -A,) (A.4) 

we find 

I3N(N-1)+N 



X 



Np\a\ 

[ d/xi • ■ • d/ijv I A(/z) | 2 ^ e" E ^ w e" 31 ^ E£l m? 

(A.5) 
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When a — > — oo, the integral is of order 1 and converges to a Selberg integral with 
Laguerre-type weight: 



d// 1 ---d^|A(//)| 2 ^e-^ w 



(n£ir(i+jfl 
r(i + /3) iV r(i + iV/3) 



(A.6) 



Thus, we have: 



lim 

o— ¥— OO 



hx(X Nt p(a))-[0N 



2t 



f3N 2 — + hit + {p - l)N (1 - ln(27r) - hit) 



ln]a| [• + ((3 - l)JVln|a| 



-(3N 2 lnN + (3N 2 ( --ln/3 ) + (/3 - 1)N In JV + (/3 - 1)N (-1 + ln(2vr) + In 0) 



N 



-lnT(l + N(3) - iVlnr(l + 0) + ^21nT(l + 



(A.7) 



Hence: 

Cat,/? 



/3A^ 2 lnA^ + /3A^ 2 - -ln/3 - l)jVlniV 



exp 



5 ^ (nf =1 r(i+j/3))' 

+(/? -l)AA(-l + ln(27r) + In • 21-^1- V 



r(l + 0)*fY{\ + AT/3) 

(A.8) 

Note that this normalization constant is the correct one when there is one hard edge. 
The result is different when there is no hard edge. 

B Reminder on special functions 
B.l Gamma function 

We recall the definition of the Gamma function. For Res > 1: 



2—1 —a 



T(z) = I daa z ~ L e 
'o 



(B.l) 



It satisfies the functional equation: 



T{z + l) = zT{z) 



(B.2) 
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It can be analytically continued on the complex plane, and 1 /T(z) is an entire function, 
which can be defined by the product: 



1 oo 

r(z) ' m 

v ' m=l 



J] 1 + - )e-* (B.3) 



Its asymptotic expansion is given by Stirling formula: 



, ,-,/ % -, In z ln(27r) , , JN 
In T(z + 1) = z In z - z + + — ^ — i + o(l) (B.4) 



It satisfies the reflection property: 



r(i - *)r(z) = — ^ (b.5) 



and the "addition of angles" formula: 

p-i / \ 

JJ T ( - + 2 J = (2tt) ^ p-^+s r(pz) (B.6) 

m=0 ^ ^ ' 

B.2 Riemann zeta function 

We recall the definition of the zeta function. For Res > 1: 

oo 1 

coo - E-i ( R7 ) 

It can be extended by analytical continuation on the complex plane, by the reflection 
relation: 

r (0 vr-/ 2 coo = r (^) 7T- V c( i _ s) (B . 8 ) 

B.3 Barnes G-function 

We recall that the Barnes G-function [4] is an entire function satisfying: 

G{z + 1) = T{z)G{z) (B.9) 
At integers, it coincides with the product of factorials: 

fc-2 

VfceN G(k) = (B.10) 

We take the convention that a void product is equal to 1. Its asymptotic expansion for 
z — > oo is a classical result [27]: 

, -\ \ o /In 2; 3\ ln(27r)2; h\z 

lnG(z + 1) = z 2 ^— - -J + - — + C'(-l) + o(l) (B.ll) 
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We quote some special values and relations [1]: 

lnG(1/2) = + ^ (B, 2 ) 

lnG(1/ 4) + ln G (3/4) .^_MM,!^/l) + £«zl) (B , 3) 

, , s , , s 71n3 ln(2vr) lnr(l/3) 8C'(-l) 

In G 1/3 +lnG2/3 = — { —t ^ + B.14 

3d 3 3 3 

ln G (l/6) + ln G (5/6) = _^_^ + !^_li^ + 5CH) 

(B.15) 

In principle, G(p/q) can be expressed in terms of an increasing number of fundamental 
constants when p and q increase. We do not try in this article to express G(p/q) in 
other ways. We mention an "addition of angles" formula [17]: 

X> G g + z + i) = ( p -I) C ' ( - 1) + i!l^±i)- zln r W 



m=0 



9 1 \ lap ( vn \ , „ /m 



(B.16) 



C Large AT expansion of some Selberg integrals 

We have defined: 



/ dA 1 ...dA^|A(A)|^e-^S- A - (C.l) 
Jr n 



This is a Selberg integral, and it is given by: 

f \ 2 2 IL=i r ( 1 +.7/ 3 ) 



= (2vr)^ 2 



N/3 J T(l + /3) 



A 1 



In the main text, we are interested in: 



(C.2) 



V{a) = = ^ eFia) ( C - 3 ) 

up to a factor equivalent to 1. So, we need the large N asymptotics of lnX/v^ and 
lnCjv,^ up to o(l). 
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C.l Asymptotics of J2f=i In T(l + j/3) 



j-- 



We first determine the large N asymptotics of Y2f=i l n T(l + jP) for a general (3. We 
start with an integral representation of lnT: 

\nT(z) = ( z -l-)lnz - z + h 



We have: 



2 7 2 



+ / ^— — -i + tt) (C4) 



iv 



^lnr(l+j/3) 

7=1 
iv 

^ln(j/3) + lnr(j/3) 
7=1 

N 

]Tlnj + iVln(/3) 
i=i 

+ £(^)l»07^^ + A' 



ln(27r) 



7=1 



A' 

N , JV 

/3^jln(j) + -X]ln(j) + 



(C.5) 



We recall the classical asymptotic expansions 

N 



^lnj = mr(iV + l) = N\nN -N + l — + l -^l + o(l) 



7=1 

A? 



, iV 2 lniV iV 2 iVlniV IniV 1 . . . . . s 
J ^ J = — = — + —— + -- + --('(-!) + o(l) (C.6) 



2 4 2 12 12 

7=1 
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Let us determine the large N asymptotic of the integral: 

A f°° fa a 
K = V / da — — 1 + - 

1 _ _ _ r°° 



+£/. ^Vl^r- 1 



a 2 \e CT -l 2 12 

1 + - 



o 

-DO J 



_L vl 

12 Pj^j + Jo (J 2 - I - I ~'2 12 
In AT | 7B + f°°dff 1 / q i | ff 2 



12/3 12/3 7 a 2 e^-lVe CT -l 2 12, 

where 7_e is the Euler-Mascheroni constant. Collecting the results, we find: 

N 

]Tinr(i+j73) 

/3iV 2 lniV / 3 hi(/3)\ 2 /3 + 1 ri Ar 

= h B[ 1 — N 2 + - N In TV 

2 f \ 4 2 J 2 

1 /31n(/3) /3 ln(/3) ln(27r)' 

2 2 2 2 

+ /3 + f 2 1 + 3 lniV + ^ + o(l) 

where the finite part k@ is given by: 

ln(27r) /l .,, A _ 7j5 



«/3 



^(n-c'(-D) 



4 ^ V 12 7 12 /3 

00 da 1_ / a 

a 2 " e^ CT - 1 Ve CT - 1 ~ + 2 ~ 12, 



or equivalently 



m(27r) fl ... _ 



7ff 
12/3 



+ 



da 



6acoth(a/2) - 12 - a 
12a 2 (e^ - 1) 
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C.2 Asymptotics of 2jv,/3 

If we apply Eqn C.8 to the definition oHn,p (Eqn. C.2), we find: 

+ (pin(P) ~ lnr (! + ln (2vr) - lntj iV 

+ / ?+r i +3 lniV+ ^ +o(1) 

(C.ll) 

C.3 Asymptotic of Cn,p 

If we apply Eqn. C.8 to the expression of C^^ (Eqn. A. 8), we find: 
InCW,/? = (3N^N + N (-(3 + (3\n(3 + /31n(27r) -lnT(l + (3)) + f - + - ln(JV) 



6 7 



. 5-2(^ + r -) N n2 _ln^)_lng + 2Kg + o(1) 



(3 



(C.12) 



Finally we get the normalization constant for the cumulative distribution of A max . 
We have indeed 

V{a) = = ?^ eF(a) = ex P i F ( a ) + lnC V/? - InX^} (C13) 



with 



lnC„-lnI„ = ,3iV 2 (5 - ^ + N(P - 1) (-1 + ln(2x) + ^ 



+ to ^p + ^'- 3 ) + ( 5 - 2 '7^)ln2 

_MM _ m + Ks (ai4) 



C.4 Another computation of ^ for (5 rational 
C.4.1 /3 integer 



totic of Yl!j=i In T(l +j/3) when (3 is a positive integer. Here, we only focus on the finite 
term of this expansion, which we called previously Kp. Of course, one could check with 
this method that one does obtain the divergent term of Eqn. C.8. 



By definition we have Kp = ^ J=1 lnr(l + j/3) 



finite 



. Let us come back to the asymp- 



5N 



For — 1, this series is the Barnes function itself: 

N 

^lnr(l+j) = lnG(iV + 2) 

3=1 



(C.15) 



Using the asymptotic expansion of G(z + 1) given in Eqn. B.ll, we obtain for any w 
in the large z limit: 

w ln(27r) 



[]nG(z + w + l)\ 



finite 



+ ?(-!) + o(l) 



Thus, we have: 



ln(27r) 

K l = n h C (-1) 



(C.16) 



(C.17) 



As a matter of fact, one may use only the asymptotic expansion of Barnes function 
(Eqn. B.ll) to work out the case /3 integer. We start from the "addition of angles" 
formula (Eqn. B.6) and write: 



N 



3=1 



(l-P)N 0N(N+1) N 
2 P 2 2 



g G(jv + y + i) 

m=0 G ^ + lj 



(C.18) 



Then 

N 



^(iv + i)_iv' 1M 



/3-1 r 

E 

m=0 



iV+ — + 1 



m 
J 



(C.19) 



Subsequently: 



ln(27r) 



-1 r 



E 

m=0 



mln(27r) 



C'(-l)-lnG 

p-i 



m 



(^±^j ln(2vr) + /3C'(-1) - f>G (l + ^ (C.20) 



For example: 



k 2 = ^^ + 2C / (-l)-lnG(3/2) 
51n2 ln(2vr) + C'(-l) 



24 2 2 

where we used the value of G(l/2) given in Eqn. B.12 and T(l/2) = y/ix. 



(C.21) 
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Actually, Kp can be expressed only in terms of values of T at rational points, by 
means of the formula Eqn. B.16: 

| ta0 ( 1+ =).(,-i), ( _ 1)+ ^ + £= hr (=) (c, 2 , 



Thus, we have: 

Kp 

For example: 



^ + 1 W^-^-gf,rg) (c, 3) 



uina in^ ctzi) + mr(i/3) 

3 36 3 3 3 V 7 

71n2 ln(27r) C'(~l) , lnr(l/4) 
= ^- + ^ + ^ + ^^ ( C - 25 ) 

C.4.2 Value of for rational [3 

The case = p/q G Q can be treated similarly (we assume gcd(p, q) = 1). It is 
convenient to do the computation for a given congruence of N mod q. Let us write: 

N = qL + R, < R<q-1 

We have according to the duplication formula: 

N / \ p-l N , . 

n_ / .V\ /n s V-P)N pN(N+l) N -I-T -I-T_ /J m 

j=l v y/ m=0j'=l vy F 

(l-p)iV P N(N+1) N 
= (27TJ 2 p 2 9 2 



g niiG(L + § + ? + 1) ■ nu + i g (l + j + ^ 



m=0 



r I m 
9 P 



(C.26) 



Then: 



r>M)^) w ]H9 

We need the large z asymptotic expansion of In (7(1 + w + z/q) for fixed w. Starting 
from Eqn. B.ll, we obtain: 



In G ( 1 + w + - 

Q 



finite 



i-^'^ + f'M (C28) 
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Then, we have: 
ln(27r) 



K p/q 



r m 



£X> G U + 7)+^'(-i) 

m=0 r=l ^ V 



P-1 f R 

EE 

m=0 I. r=l 



1 1 

12 ~ 2 



+ E 



=R+1 



1 1 

12 ~ 2 



-i? + r m 
q p 



-R + r 



q 



m 

1 + — 

V 



+E E 

m=0 I. r=l 



-R + r m\ \ ~* 
+ — + > 

y ^ / r=R+l 



-R + r 



3-p/q-q/p p-q + 2 

12 ng " 4 



In g 



i + m^| ln(27r) 

p-1 q 



ln(27r) + mC'(-I) " E E ln G (~ + -) 

(C.29) 



Notice that the result does not depend on the congruence of N modulo q. For q = 1, 
we recover the result for (3 = p integer (see C.20). 

We give below some values of Kp. To complete the table, we have used repeatedly 
the reflection formula and the "duplication of angles" formula of the T function, and 
Eqns. B.12-B.15. 



Value of (3 




1 

2 
3 
4 


ln( f ) +C / ( i) 

5 In 2 | ln(27r) . C'(-l) 

24 2 2 
Illn3 | ln(27r) . . l n r(l/3) 

36 3 3 3 
71n2 , ln(27r) C'("l) i lnr(l/4) 
12 1 4 1 4 1 2 


1/2 
1/3 
1/4 


In 2 ! ln(27r) { C'(-l) 

2tn3 | ln(2vr) InY(l/3) C'(-l) 

5 In 2 , ln(27r) . C'(-l) , lnr(l 3 /4) 
8 1 4 1 4 1 2 


2/3 
3/2 


13 In 2 51n3 , 21n(27r) lnr(l/3) C'(-l) 

72 18 1 3 3 1 6 
ln2 | llln3 , 21n(27r) lnr(l/3) C'(-l) 
4 72 1 3 3 1 6 
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